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Abstract 

In this note, we derive (to third order in derivatives of the fluid velocity) a 2+1 
dimensional theory of fluid dynamics that governs the evolution of generic long- 
wavelength perturbations of a black brane or large black hole in four-dimensional 
gravity with negative cosmological constant, applying a systematic procedure de- 
veloped recently by Bhattacharyya, Hubeny, Minwalla, and Rangamani. In the 
regime of validity of the fluid-dynamical description, the black-brane evolution 
will generically correspond to a turbulent flow. Turbulence in 2+1 dimensions 
has been well studied analytically, numerically, experimentally, and observation- 
ally as it provides a first approximation to the large scale dynamics of planetary 
atmospheres. These studies reveal dramatic differences between fluid flows in 
2+1 and 3+1 dimensions, suggesting that the dynamics of perturbed four and 
five dimensional large AdS black holes may be qualitatively different. However, 
further investigation is required to understand whether these qualitative differ- 
ences exist in the regime of fluid dynamics relevant to black hole dynamics. 



1 Introduction 



A particularly fascinating aspect of gauge theory / gravity duality is that the high- 
temperature deconfined phases of gauge theories are mapped via the correspondence 
to black hole or black brane geometries on the gravity side pQ. This enables analytic 
calculations of equilibrium properties of strongly-coupled gauge theory plasmas, such 
as the deconfinement temperature or the high-temperature equation of state, via a 
classical analysis of the corresponding geometries. 

In the past several years, the connection between gauge theory plasmas and black 
geometries has also been exploited in the near-equilibrium regime. On the gauge theory 
side, it is known that long-wavelength fluctuations about the equilibrium state at high 
temperatures should be described effectively by fluid dynamics equations. Specifically, 
near-equilibrium configurations may be characterized by a local temperature and a 
local fluid velocity w M , in terms of which the stress tensor of the theory can be written 
as a derivative expansion, 



The fluid dynamics equations are then simply the local conservation equations for this 
stress tensor 



The information that distinguishes different field theories in this regime is the set 
of coefficients appearing in the derivative expansion of the stress tensor, namely the 
energy density and pressure at zeroth order, the shear viscosity and bulk viscosity at 
first order, etc... . 

While the energy density and pressure are equilibrium quantities, which may be 
read off from the dual Euclidean black-brane geometry, the viscosities and all higher 
order coefficients are only relevant for time-dependent processes, and their computation 
requires a real-time analysis. The first calculation of viscosity appeared in [3J, where 
the authors were able to determine the viscosity of maximally supersymmetric gauge 
theory plasma using a Kubo-type formula, i.e. by studying the low-frequency limit of 
real-time correlation functions for the stress tensor. Following this, there have been 
a host of works exploring the fluid dynamics description of high temperature gauge 
theories with gravity duals (see p2] for a recent review), motivated in part by the fact 
that the high-temperature behavior of these theories appears to be qualitatively similar 
to real QCD. Indeed, this AdS/CFT analysis currently provides the best understanding 
of certain properties of the quark-gluon plasma produced experimentally in heavy-ion 
collisions at RHIC. 

Recently, Bhattacharyya, Hubeny, Minwalla, and Rangamani [2] have developed 
a more direct approach for determining the effective fluid dynamical description of a 
near equilibrium gauge theory plasma with a weakly curved gravity description. This 
approach provides a systematic way to determine the coefficients in the derivative 
expansion of the stress-energy tensor to an arbitrary order, and the explicit gravity 
solution dual to a given solution of the fluid dynamics equations. 



= p(T)rT + (e(T) + p(T))uV + 0(du, &T) . 





(2) 
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The idea of [2] is simply to start with the black brane solution dual to an equilibrium 
plasma with uniform temperature T and velocity [3, allow the temperature and velocity 
parameters in the metric to be arbitrary slowly-varying functions of the field theory 
directions, and try to add corrections to this metric order by order in the derivatives of 
temperature and velocity so that the corrected metric is a valid solution to Einstein's 
equations. This perturbative procedure works subject to certain constraints on the 
temperature and velocity fields; these constraints turn out to be exactly the conser- 
vation relations ([2]) for a particular choice of coefficients in the derivative expansion 
(CQ) of the stress-energy tensor. In [2J this procedure was carried out in the case of 
five-dimensional gravity to determine the complete set of fluid dynamics coefficients 
up to second order in the derivative expansion. The result for the viscosity matched 
previous computations, and a subset of the second order coefficients they compute are 
reproduced using correlation function methods in |4i], which appeared simultaneously. 

1.1 Results 

It is worth emphasizing that the calculations in [2] are purely within the context of 
Einstein gravity with a negative cosmological constant, and do not assume in any way 
the correctness of gauge theory / gravity duality or any connection to a dual field theory. 
Thus, the results are also interesting from the perspective of better understanding 
classical gravitational physics. To reiterate, they say that generic long-wavelength 
perturbations around a black brane solution in AdS (or, as we will discuss below, a large 
AdS black hole) are described by fluid dynamics equations. Of course, nonlinear fluid 
dynamics is extremely difficult to study in general, but decades of analytical, numerical, 
and experimental research have provided us with a good qualitative understanding and 
some general quantitative laws for the behavior of fluids in various regimes. We can 
therefore hope to apply what is known about fluid dynamics to learn new qualitative 
and quantitative lessons about the generic behavior of gravity in this nonlinear regime. 

To apply these ideas to the most familiar case of 3+1 dimensional gravity, a first 
step is to generalize the work of [2 J (which focused on 3+1 fluid dynamics and therefore 
4+1 dimensional gravity) to one lower dimension. This is the main technical goal of 
the present paper. 

Before giving our results, it will be useful to define exactly what we mean by the 
temperature field and the velocity field associated with a given field theory (or more 
generally, with a given stress tensor). For a given point x in the fluid, there will be 
some velocity /3 l (x), such that an observer at x moving with this velocity will observe 
T 0l (x) = 0. We define the proper velocity for this point to be u^(x) = (1, /3(x))/(l — 
/3 2 (x)) 1//2 , and we define the temperature for this point to be the temperature observed 
by the comoving observer (the temperature related by the equilibrium equation of state 
to the energy density T 00 (x) as measured in the comoving frame). This definition of 
velocity may be expressed covariantly as 

u„(T^ + u v u a T m ) = 

or equivalently u^T^" = where T represents all terms in T involving derivatives. 
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With these definitions of velocity and temperature, our result is that the evolution of 
long- wavelength perturbations to the uniform black- brane solution of 3+1 dimensional 
gravity with negative cosmological constant is described by the conservation equations 
(j2J) where the stress energy tensor to second order in the derivative expansion is 



1 /4vrT\ 3 , „„ „ „ „ s /4ttT x 2 



= g ( — ) W + 3uV) - ( — J a*™ (3) 



1 /4vrT 
+ 18 



(V3tt - 9 ln(3) + 18)Ef + 2(^71 - 9 ln(3))Ef 



The terms beyond leading order are traceless symmetric tensors which vanish when 
contracted with u^. Defining the operator 

= rf v + u'V (4) 

which projects vectors in the directions orthogonal to u and 

= \pSPp + \p$K - I^Pop ( 5 ) 

which project tensors into traceless symmetric tensors orthogonal to u^, we have 

o*> = n^(sV) 

uf" = -P%Pp {d a u p - d p u a ) 

= K P {^ X ) (6) 

where D = u a d a . From the expression above, it is straightforward to check that the 
ratio of viscosity (—1/2 times the coefficient of a) to entropy density (which works out 
to 1/T times the coefficient of u^ii") comes out to the universal [HI [T3, [TB] value of 
l/(47r), first shown in this case by [5]. 

As reviewed in (HE], the AdS isometries (or equivalently, the conformal symmetry of 
the dual fluid) require that the fluid stress tensor is traceless and covariant under a Weyl 
transformation of the boundary metric (i.e. the metric of the space on which the dual 
fluid lives, which we take to be 2+1 dimensional Minkowski space). The two structures 
Si and E 2 turn out to be the only traceless u-orthogonal Weyl- covariant structures at 
two-derivative order relevant for a flat-space theory. In higher dimensions, there are 
four such structures [H [8], but two of these (the ones associated with parameters Ai 
and A3 in the notation of [1]) vanish identically for 2+1 dimensional fluid dynamics. 
In the notation of [3], our results for the second order coefficients are 

Tn = ^(^-9111(3) + 18) 
A 2 = — (V37r-91n(3)) 
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The result for r n agrees exactly with the calculations of [6j [7] , who computed this using 
real-time correlation functions, while our result for A2 is new. 

In addition to our expression for the stress tensor, which through the conservation 
equations gives the equations for the effective fluid dynamics to third order in the 
derivative expansion, the other main technical result is the equation (1211) . which gives 
the explicit metric dual to an arbitrary solution of the fluid dynamics equations, valid 
to second order in derivatives. 

1.2 Possible qualitative differences between four and five di- 
mensional black hole evolution. 

The derivative expansion we have presented should be a good description when the 
length scale L of perturbations is much larger than the inverse temperature 

LT > 1 . (7) 

In this limit, the most important terms are the leading terms in the derivative ex- 
pansion, and these are the terms conventionally included in studies of fluid dynam- 
ics. While our calculations explicitly refer to perturbations around a uniform black- 
brane solution, dual to an unbounded fluid, the same fluid dynamics on a sphere of 
radius R ^ should describe the evolution of perturbations around a large AdS- 
Schwarzschild black hole (i.e. one whose radius is large compared with the length scale 
of AdS curvature). 

At the level of the equations of motion, our results look very similar to the ones 
obtained in [2] for 4+1 dimensional gravity. Despite the similarity of the equations, the 
relation to fluid dynamics suggests that there could be dramatic differences between the 
physics of large AdS black holes in four versus five dimensions. To understand these, 
we first note that the regime LT ^> 1 where our derivative expansion is valid is also 
the regime of large Reynolds number (where the viscous terms are small relative to the 
leading order terms), so we expect that generic evolution of the long- wavelength per- 
turbations about the black brane will correspond to fluid flows in the turbulent regime. 
In 3+1 dimensions, turbulent flows are characterized by an "energy cascade" in which 
large scale eddies give rise to smaller scale eddies, tending towards configurations of 
increasing disorder and transferring energy down to scales where viscosity becomes im- 
portant and energy is dissipated (converted to heat). Thus, 3+1 dimensional fluids are 
rather efficient at dissipating initial fluctuations, regardless of how small the viscosity 
is. For a basic review of fluid dynamics and turbulence, see P, \TU\ . 

In contrast, turbulent flows in 2+1 dimensional fluids are characterized by an "in- 
verse cascade," in which smaller scale eddies merge into large scale eddies, eventually 

1 This follows from a simple AdS/CFT argument, since a thermal CFT on Minkowski space dual 
to the black-brane arises from the same CFT at finite temperature on a sphere in the limit of large 
sphere volume and fixed temperature. The CFT on sphere with R 3> 1/T is dual to a large AdS black 
hole, and this CFT should have locally the same long-wavelength physics as the flat-space one, since 
correlation lengths in this limit will be much smaller than the sphere size. 
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creating large persistent vortical structures from which the energy is dissipated only 
very slowly, since the dissipative terms are relatively unimportant at large scales (see 
section 9.7 of [10], and also [H]). In unbounded systems the late-time behavior for 
generic initial conditions is believed to be a "dilute gas of vortices" while in bounded 
systems (such as fluid flow on a sphere), the flow evolves to a single vortex or a small 
number of vortices with size comparable to that of the system. These two-dimensional 
phenomena are relevant in nature, since two-dimensional turbulence provides a first ap- 
proximation to large scale motions in planetary atmospheres, and also to the evolution 
of large-scale oceanic currents (examples of the persistent vortical structures include 
intense tropical storms, the polar vortices, and the Great Red Spot of Jupiter). 

These qualitative differences between two and three dimensional fluid dynamics 
suggest that the evolution of generically perturbed AdS black holes may exhibit signif- 
icant qualitative differences between two and three dimensions. Taken at face value, the 
results of the previous paragraphs would suggest that black brane and large black-hole 
configurations in five- dimensional gravity with negative cosmological constant equili- 
brate much faster than those in four dimensions, and furthermore, that the late-time 
behavior for generic perturbations around a four- dimensional black brane or large AdS 4 
black hole will be characterized by the formation of long-lived large-scale vortical struc- 
tures, related to the more familiar ones observed in planetary atmospheres . 

However, there is an important caveat U The fluid dynamical results we refer to 
are valid in the context of Navier-Stokes equations for non-relativistic fluids which 
have p <C e (in analytical and numerical studies, the fluids are usually also taken 
to be incompressible). On the other hand, the fluids we discuss have pressure and 
energy density of the same order of magnitude; in other words, they are microscopically 
relativistic. For such fluids, the equations of motion differ from Navier-Stokes equations 
even in the limit where all the macroscopic velocities are small. It is plausible that the 
qualitative differences between two and three dimensional fluid dynamics are not just 
restricted to the Navier-Stokes regime. However, before making any specific conclusions 
about about differences between four and five dimensional black hole evolution, it will 
be important to understand whether the qualitative differences we have noted (or other 
differences) exist in the regime where the fluids are microscopically relativistic. A recent 
discussion highlighting some features of two-dimensional turbulence in relativistic fluids 
has appeared in [T7] . 

It is interesting to ask whether any of this discussion is applicable to astrophysical 
black holes (which should behave like small black holes in AdS). Unfortunately, the 
present analysis doesn't shed any light on this question, since for ordinary Schwarzschild 
black holes (or small AdS black holes), the length scale associated with temperature 
is just the size of the black hole. Thus, the condition (J7]) for the validity of the hydro- 
dynamic approximation cannot be satisfied; equivalently, all higher order terms in the 
expansion of the stress-energy tensor are equally relevant. 

2 I would like to thank Veronika Hubeny and Mukund Rangamani for motivating me to understand 
this better. 
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1.3 String theory applications 



While we have emphasized that our results are apply completely within the context of 
four- dimensional gravity, there are also important applications within the context of 
string theory using the AdS /CFT correspondence. According to the correspondence, 
the fluid dynamical theory we have worked out provides the effective description of 
long-wavelength perturbations around the thermal equilibrium state for any 2+1 di- 
mensional conformal field theory with a weakly curved gravity dual, provided that 
Einstein gravity with negative cosmological constant is a consistent truncation of the 
corresponding gravity theory. The simplest example is the SO (8) conformal field the- 
ory that provides the worldvolume theory of a large number of M2-branes in M-theory. 
For this theory (and other theories) there are locally conserved charges in addition 
to the energy and momentum. The fluid description we have found applies to neutral 
fluctuations about the uncharged equilibrium state, but it should be straightforward to 
determine a more complete effective theory which allows for nonzero charge densities 
by applying a similar analysis starting with a charged black-brane solution. 



The remainder of this paper is organized as follows. In section 2, we present the 
perturbative calculation of the metric for a generic non-singular perturbation of the 
3+1 dimensional black-brane solution to second order in the derivative expansion, and 
derive the associated constraints on velocity and temperature fields. In section 3, we 
use our expression for the metric to calculate the boundary stress tensor, using the 
prescription of [T3]. This stress tensor is interpreted as the stress tensor of the dual 
fluid, and we find that the constraint equations found in section 2 coincide with the 
conservation equations for this stress tensor, expanded to second order in derivatives. 

2 Construction of the perturbed solution 

In this section, we apply the methods of [2] to perturbatively construct a solution of 
Einstein gravity with negative cosmological constant dual to an arbitrary solution of 
particular 2+1 dimensional fluid dynamics equations that we will also construct in the 
process. 

2.1 Unperturbed metric 

We begin with Einstein's equations in 3+1 dimensions with negative cosmological con- 
stant 



1.4 Outline 



Ru 



- -guR - Xgu = , 
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where for convenience, we choose units such that A = 3. The uniform black brane 
solution corresponding to a temperature T is^| 

ds 2 = 2dvdr — r 2 f(br)dv 2 + r 2 dx 2 

where3 

Boosting this solution to a proper velocity along the x directions, we obtain the 
three-parameter family of solutions 

ds 2 = -lu^dr - r 2 f{br) Ufl u u dx^dx u + r 2 P iiV dx lL dx v (9) 

where is the proper velocity 



l-P \/l-/3 2 

and Pfj, u = r]^ + Ufj,u u is the projection into the directions orthogonal to u, defined so 
that u^P^ = 0. 

This solution depends on the temperature T and velocity (3, so we have a map 
between equilibrium configurations of a fluid moving at constant velocity and uni- 
form black-brane solutions. We would now like to extend this correspondence to near- 
equilibrium configurations in which the temperature and velocity of the fluid are allowed 
to vary in space, evolving in time according to fluid dynamics equations. Thus, we now 
promote b and (3 to general functions of space and time (assumed to vary slowly on the 
scale 1/T) in the metric (jUJ), and try to add corrections to this, order by order in deriva- 
tives of b and (3, such that the corrected metric is a solution to Einstein's equations 
(jSJ). As in we will find that this is possible as long as b(x^) and /3(x^) satisfy fluid 
dynamics equations with particular choices for the pressure, energy density, viscosity, 
and higher order coefficients. 



2.2 Overview of perturbation theory 

We now review the perturbative procedure developed in [2]. At each order, we begin 
with a metric g( ra_1 ) depending on b(x^) and which solves Einstein's equations 

up to terms involving n — 1st derivatives of b and /3, but for which terms in Einstein's 
equations involving nth and higher derivatives do not cancel. For n — 1, the metric is 
just the one in (jH]) with b and /3 taken to depend on space and time. We then try to 

3 Thcse coordinates, introduced in [2] are related to more typical coordinates where the metric is 
diagonal by a coordinate transformation v = t + h(r), where h'(r) = l/(r 2 /(r)) 

4 As usual, the relationship between temperature and the parameter b may be determined by 
demanding that there is no conical singularity in the Euclidean continuation when the Euclidean time 
direction is chosen to have periodicity 1/T. 
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add an nth order metric g^ 1 ' such that the complete metric solves Einstein's equations 
up to nth order in derivatives. 

Throughout the perturbative procedure, we work with the gauge choice [2] 

g rr = g r , oc u, (<? (0) r <£ >0) = 2gff + ±g<p = (10) 

The metric is further constrained by the requirement that its components g^ r and 
g^u must transform as a vector field and a tensor field under 2+1 dimensional Lorentz 
transformations where and b are taken to transform as vector and scalar fields. This 
follows from the covariance of the leading order metric and the covariance of Einstein's 
equations. The most general covariant metric at nth order in derivatives satisfying our 
gauge condition is 

(ds 2 Y n) = \u ll u v dx ll dx v -2h n u^dx fl dr - r 2 h n P ViV dx^dx v 
T 2 

--{j^vUpdxrda? + r 2 (a n )^dx^dx u (11) 

where we may choose 

u„f = UypT = -q^oT = . 

Here k n (r,u tl (x),b(x)), h n (r,u ll (x),b(x)), j^(r,u IM (x),b(x)) and a^ / (r,u IJ ,(x),b(x)) are 
expressions involving a total of n derivatives of u^(x) and b(x), that transform respec- 
tively as two scalar fields, a vector field, and a symmetric traceless tensor field when 

and b are taken to transform as vector and scalar fields under 2+1 Lorentz trans- 
formations. Note that the only dependence of these functions on the coordinates x M is 
through the functions 6(x M ) and u(x^). 

In order to determine the functional dependence of k, h, j, and a on r, b, and u, it 
is enough to solve Einstein's equations at any particular point, so we choose to work 
at x^ = and further, take coordinates where 6(0) = 1 and /3j(0) = 0. 

In this case, the extra terms ( ITTi) that we add at nth order in perturbation theory 
reduce to 

(ds 2 ) in) = ^p-dv 2 + 2h n {r)dvdr - r 2 h n {r)dx l dx l + -ji^dvdx 1 + r 2 a^(r)dx i dx j 
Using this undetermined expression and the results from lower order, we now take 

g = g^ + ...g^\ 

Taylor expand about x = to nth order in derivatives of b and f5, plug into Einstein's 
equations, which may be rewritten as 

W u = R u + 30/ j = , (12) 

and try to choose h n , k n , j l n) and a l J so that all terms with n derivatives cancel (terms 
with less than n derivatives will cancel assuming that we have correctly carried out 
the perturbative procedure at lower orders). We will find that this is possible so long 
as b(x) and (3(x) obey certain constraint equations, which amount to the equations of 
fluid dynamics with specific coefficients for pressure, energy density, viscosity, etc... 



S 



2.3 General structure of perturbation theory 

Carrying out the procedure above, we find that at nth order in perturbation theory, 

and ajf may be determined by the four equations 



the functions h n , k n , j 



r 4 /(r)PK ( r n) - W % 



rr 
(n) 
ii 



w 



W; 



(n) 



(n) 
H 










dr 1 



— k r . 



-( 



+ (l-4r 3 K(r))) 



in). 



(»)/ 



Hi 



S3 



where Wj™ denotes all n derivative terms in Wjj (defined in f fT2l) ) and the equations 
on the right come from dividing these terms into those that arise from g( n > and those 
that arise from lower order terms in the metric. 

The general solution to these equations is given by the specific solution 



K{r) 
k n {r) 
jn(r) 



% Hdyy^iy) 
dxSt\x) + r -{l 



oo poo 
2 



J x 



dxx'J jSr(y) 



(14) 



plus the general solution to the homogeneous (source-free) equations, 



K{r) = s n + 



k n (r) = u n r + — -2rs n 

3n{r) = a n + —r 

I i 

«n = c n + -d n ln(r 2 ) . 

6 r 



(15) 



With our choice for the the specific solution, it turns out that we can/must set all 
coefficients in the homogeneous solution to zero. First, we must set s n , b n , c n and d n 
to zero in order to get a nonsingular solution that preserves the original asymptotics. 
The constant t n may be set to zero by a coordinate transformation r — > r + t n / (2r 2 ) 
that preserves our gauge choice. The coefficients u n and a n may be adjusted freely by 
redefinitions of b and (3. But we chose to fix the definitions of b and (3 by demanding 
that u l j L T llv = 0. For any non-zero u n or a n , this condition is violated, so we must 
demand that these constants vanish. Then the b and (3 that appear in the metric are 
directly related to the inverse temperature and velocity of the fluid. 

Once we have solved for h, k, j, and a at a given order, the remaining four equations 
W rv = W vv = W V i = give rise to constraints on the set of n-derivative expressions 
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built from temperature and velocity. These constraints turn out to be precisely the 
fluid dynamics equations 

Taylor expanded to nth order in derivatives about the point x — 0. The constraints 
at order n involve the stress tensor at order (n — 1), which can be computed as the 
boundary stress tensor of the metric obtained at the previous order in perturbation 
theory. 

Given the solutions for h, k, j, and a we can write the metric at order n in the 
derivative expansion. From this, we can calculate the boundary stress-energy tensor 
(which we interpret as the stress-energy tensor of the dual fluid) at order n, using the 
prescription of [13] reviewed in section 3. 

In the next two subsections, we provide explicit details of the perturbative calcula- 
tion at first order and second order in derivatives. 



2.4 First order 

At first order, the source terms in equation ([TBI are 



s« = o 

1 



s; (1) = -~d v p 

S<?> = 2ra 

where ^ ^ 

a ij = -jidifii + djA - -SijdkPk) . 

With these sources, the metric components at first order may be calculated using 
and give 

hi(r) = 
h(r) = r 3 di(3i 
Ji(r) = r 2 d v /3 



2 2 



aj = 2aF(r) ~ a(- - — + 0{r~ A )) 
r 6r 6 



where 




Fir) = -^Tan" 1 [ ^(2r + 1) ) + ^ ln(l + i + 1) + ^ . (16) 



With these assignments, we find that the remaining Einstein equations are satisfied if 
and only if 

d v b = ^difr 

d l b = d v (3 (17) 
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These are the terms obtained at first order in derivatives in the conservation equation 
dfjT^" = 0, where we need only keep the non-derivative terms in the expression ([3]) for 
T. 



2.5 Second order 

To discuss the results at second order, it is useful as in [2] to catalogue all the scalar, 
vector, and traceless symmetric tensor terms built from b and (3 at two-derivative order, 
since these are the expressions that will appear in the second order metric. We find: 



Scalar 


Vector 


Tensor 




Si = 


dlb 


V u 


= d v dib 


T U j 


= d^b - \8 i: 


,d 2 b 


s 2 = 


dfb 


v 2l 


= d 2 A 


%ij 


dy&ij 




5 3 = 


d v diPi 


v 3t 


= didjPj 


%ij 


dvftid v j3j 


\8ij(d v j3k) 2 


s 4 = 


d v Pi d v Pi 


v 44 


= d 2 (3 l 


%ij 


= dkPidkPj — 




s 5 = 


(<W 


v 54 


= dvpidjpj 


T 


= dif3kdjf3k — 


\5 l3 (d k fr) 2 


s & = 




v 64 


= dvPjdiPj 












v 74 


= dvPjdjPi 









Here, we have excluded pseudoscalar, pseudovector and pseudotensor terms, since these 
will not appear in the metric. Also, note that in 2+1 dimensions, we have 

1 1 

dk(3(idj)Pk — ^ijdkPidiPk = + 

so we do not need to include the structures on the left, which are independent in higher 
dimensions. 

The Einstein equations at second order in derivatives have a solution provided that 
the following constraints are satisfied 

Si = ~ 2^ 4 + 4^ 5 
iS 2 = S 3 + S 4 — -Sq + 1S7 

v 2 = ^v 3 + iv 5 -^v 6 -v 7 

Tx = T 2 + % + -T A + -T b 
1 

d v diipj] = -dkPkdttPj] . (18) 



In addition, the first order constraint (fTTj) must be corrected to include higher order 
terms 

d v b = ^di/3i - ^S 7 
11 



dib = d v (3- l -V A - 2 -Vr + 2 -V 6 . (19) 



It may be checked that the four constraints (119]) and the nine constraints (fl8l) are 
exactly the 4+9 equations arising from the equations d^T^ v \ x= Q = and d a d^T^ u \ x= Q = 
expanded to second order in derivatives]^] 

Taking into account the constraints, we can now write the second order source terms 
in equation ffl3|) in terms of the independent two-derivative expressions and solve for 
the metric at second order. We have 

Sf = 2 «S3 + i<S5-2(l + i)<S 6 + F 2 (r)«S 7 

S j = 2^ Vs ~ 2r 2 (l + r + r 2 ) V ' ~ 2r 2 (r 2 + r + 1) ^ 6 ~ V ^ 

= F 3 (r)(T 2 + T 3 ) + F A (r)% + F 5 (r)T 5 (20) 



where 



F l(r) = W + *l, F<r) 



, s 1 + 3r + 4r 2 - 4r 3 - 6r 4 - 4r 5 4r 3 + 4r 2 + 2r - 1 
F 2 r = 2F r ~ + t — 

F 3 (r) = 2rF(r) - ^ + 1} 
r A + r + 1 

F 4 (r) = --F(r x ' 



2" v ' 2(r 2 + r + 1) 

. , 3 , , 2r 2 + 2r - 1 

Fz\r) = —rF[r) — 

5V ; 2 K J 2(r 2 + r + 1) 

and F(r) was defined in f[TBT) . 

With these sources, the metric components at second order are given by the expres- 
sions in (1141) . For our purposes of calculating the second order stress tensor, we only 
need to know the asymptotic behavior of the various functions for large r. We find: 

r A 2 4 
fci(r) = r 2 (S 7 + ±S 6 -±S 5 -S 3 )+0(r ) 

Ji(r) = K-^V 3 - \V, + ^V 6 ) + Oir- 1 ) 

«i = ^(T5-r 4 ) + i(l(^-ln(3) + 2)(4r 2 + 4T 3 -r 4 + 3r 5 ) + i(T 4 -T 5 ) 



5 In [2], the temperature field b was decomposed into terms such that the constraints set 
term b^ n ' equal to a sum of n derivative terms involving (3 without any further corrections. Our b 
corresponds to their J2 n b^ n \ 
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+0(r~ 4 ) 



In all these expressions, the higher order terms do not give any finite contributions to 
the boundary stress tensor. 

2.6 Covariant form of the metric 

We have now derived complete expressions for the functions h, k, j, and a appearing 
in the expression ffTTj) for the metric corrections at first and second order in derivatives. 
However, to derive these expressions, we have been working at a particular point x = 
with a choice of coordinates where 6(0) = 1 and /3j(0) = 0. To recover the general 
expression for the metric without these assumptions, we only need to rescale coordi- 
nates r —>■ br, x p —>■ x^/b and find Lorentz covariant expressions k, h, j^, and a pu that 
reduce to our expressions above in the frame where /?«(0) = 0. 
The final result for the metric is 

ds 2 = —2u^dx^dr — r 2 f{br)u^Uydx^dx y -\-r 2 P^ v dx^dx u 

+rd\u x u fM u u dx /J 'dx u — ru x d\{u fJi u v )dx p dx 1 ' + 2r 2 bF (br)a py dx p dx v 



^^-u p u v dx p dx v - 2b 2 h 2 (br)u p dx p dr - rWh^P^t 



2 

— —j2u(br)u fJj dx fJ 'dx l/ + b 2 r 2 a 2 u (br)dx p dx u (21) 

where k 2 , h 2 , j 2 and a 2 are defined in terms of the sources (1201 by (TH|) . and we replace 
the expressions S, V, and T built from derivatives of (3 with the covariant ones given 
by 

5 3 = u p d p d u u u -Vu»Vv? 

5 4 = Vu p Vu p 

5 5 = {dXf 

v 3 „ = p;(w-w) 
v 4M = p; L p pv d p d v u a 

V 5/1 = P°d u u v Vu a 
V 6/1 = P;Vu u d a u u 
T£ v = U^(Vd a u ) 

Xf u = Yl^\d p u a d p u p + Vu a Vu p ) 

rf = ng(0V0%) 



Here, P pu and were defined in (jlj) and ([5]) and V = u a d a 
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3 Stress tensor 



The boundary stress-energy tensor (which we associate with the stress-energy tensor 
of the dual fluid) may be computed using the prescription of [13], as 

= lim (r 5 (0^ - 7^9 - 2 7 ^ - 

r— >oo 

Here, is the extrinsic curvature for the surface of constant r, which may be calcu- 
lated as 

e IJ = v V - «VW 

where v K is the vector field of unit vectors normal to the surface of constant r, deter- 
mined by 

gi^v 1 = o gijv 1 ^ 1 = l . 

All other tensors are constructed using the boundary metric 7^ induced on the surface 
at fixed r, and G^ v is the Einstein tensor calculated from this metric (with indices 
raised by 7). 

Using the metric we have derived, it is straightforward to calculate the stress tensor 
to second order in derivatives. Note that one derivative and two derivative terms in 
the metric g^ v of order r and 1 respectively give rise to potential divergences in the 
stress-energy tensor (scaling as r 2 and r respectively), but these all cancel. The only 
terms contributing to the finite part of the stress tensor are the 1/r terms in g^ V) which 
include the zero-derivative terms proportional to u^u u and the one and two-derivative 
terms in the symmetric traceless spatial terms in the metric (i.e. the cty terms). The 
final result for the stress tensor to two- derivative order is given as equation (J3j). It may 
be checked that the structures appearing at second order may be written in terms of 
the covariant two-derivative tensors defined in the previous section as 

Si = T 2 + % + -T b + -T 4 
1 1 

2-jo — — Ik — —1a'. . 

1 4 & 4 ' 
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